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SUMMARY 


The  lateral- directional  dynamic -response  characteristics  of  a  large 
flexible  airplane  are  presented  for  flight  conditions  of  0.6  to  0.82  Mach 
number  and  an  altitude  of  35*000  feet.  Transient  responses  at  various 
stations  on  the  airplane  to  aileron  and  rudder  pulse-type  disturbances 
have  been  recorded  and  converted  into  frequency-response  form.  A  fre¬ 
quency  range  of  0.5  to  25  radians  per  second  is  presented,  which  includes 
the  Dutch  roll  and  rolling  modes  and  also  three  structural  modes.  The 
measured  frequency  responses  are  compared  with  those  obtained  from  trans¬ 
fer  functions  based  on  the  lateral-directional  equations  of  motion,  which 
include  first-order  effects  of  structural  deflections.  Reasonably  good 
agreement  is  obtained  in  the  low  frequency  range  corresponding  to  the 
lateral-directional  modes.  However,  in  the  higher  frequency  range,  the 
effects  of  the  structural  modes  are  sufficiently  large  so  that  these  pre¬ 
dictions  become  inadequate  even  for  the  quantities  measured  near  the 
center  of  gravity.  Comparisons  are  also  made  with  frequency  responses 
obtained  from  transfer-function  coefficients  calculated  for  the  rigid 
airplane  to  indicate  the  effects  of  flexibility. 

Calculated  simplified  transfer-f unction  coefficients  based  on  the 
lateral-directional  equations  of  motion  are  compared  with  coefficients 
obtained  by  matching  transient  time  histories  with  an  analog  computer. 
Aerodynamic  derivatives  extracted  from  measured  transfer -function  coef¬ 
ficients  are  also  compared  with  estimated  derivatives.  In  general,  good 
agreement  between  measured  and  estimated  values  was  obtained. 

INTRODUCTION 


Because  of  the  increased  importance  of  structural  flexibility  on 
aircraft  stability,  control,  and  loads,  the  NACA  is  currently  flight 
testing  a  large  flexible  35°  swept-wing  airplane  over  its  entire  operating 
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range  of  flight  conditions.  One  phase  of  this  program  is  to  investigate 
characteristics  of  the  airplane  dynamic  response  to  control- surface 
motions.  Some  results  for  the  longitudinal  case  and  a  limited  amount 
of  data  for  the  lateral-directional  case  were  presented  in  references  1 
and  2.  In  a  separate  investigation ,  reference  3;  some  measured  frequency- 
response  data  were  also  reported.  The  purpose  of  the  present  report  is  to 
present  a  more  detailed  description  of  the  lateral -directional  character¬ 
istics  of  the  airplane  at  an  altitude  of  35*000  feet. 

The  dynamic-response  characteristics  can  usually  be  separated  into 
the  relatively  low  frequency  lateral-directional  modes  and  the  higher 
frequency  structural  modes.  The  lateral-directional  modes  are  of  primary 
importance  in  defining  stability  and  control  characteristics  of  the  air¬ 
plane  as  a  whole ,  while  the  structural  modes  are  of  importance  in  the 
stability  of  the  structure  by  itself  or  in  combination  with  an  automatic 
control  system  (flutter).  In  the  present  report,  the  experimental  response 
over  the  low  frequency  range  is  compared  with  the  predicted  responses  of 
the  airplane  with  pseudostatic  effects  of  flexibility  included.  Predicted 
responses  of  the  rigid  airplane  are  shown  to  indicate  the  effects  of  flexi¬ 
bility.  The  measured  dynamic  response  is  also  examined  over  the  higher 
frequency  range  to  investigate  effects  of  the  first  three  structural  modes . 

With  regard  to  measuring  and  analytical  techniques,  the  frequency 
response  is  a  convenient  way  of  expressing  the  dynamic -response  character¬ 
istics  in  a  standard  manner  which  is  independent  of  the  particular  input 
from  the  control  surface.  Transient  responses  to  pilot-applied  pulse 
inputs  through  the  rudder  or  aileron  control  systems  were  measured  and 
transformed  to  frequency-response  form  by  means  of  the  Fourier  integral 
(e.g.,  ref.  4).  While  the  frequency  response  represents  the  dynamic 
characteristics  of  the  airplane,  a  sometimes  more  useful  form  for  the 
analysis  and  synthesis  of  automatic  control  systems  is  the  transfer  func¬ 
tion.  For  the  present  investigation,  transfer-function  coefficients 
obtained  from  the  lateral-directional  three -degree -of- freedom  case  were 
evaluated  by  matching  time  histories  on  an  analog  computer  in  a  manner 
similar  to  that  described  in  reference  4.  Corresponding  predicted 
transfer-function  coefficients  were  calculated,  with  the  coefficients 
modified  to  include  flexibility  effects,  and  are  compared  with  the 
experimental  values . 

The  measured  transfer-function  coefficients  usually  must  be  extrapo¬ 
lated  to  different  flight  conditions.  Since  this  extrapolation  can  be 
done  more  easily  in  terms  of  the  aerodynamic  derivatives,  it  is  desirable 
to  reduce  the  transfer-function  coefficients  to  derivative  form  wherever 
possible.  Accordingly,  aerodynamic  derivatives  obtained  from  simplifica¬ 
tions  of  the  principal  transfer- function  coefficients  are  compared 
with  estimated  derivatives. 

The  data  used  in  this  report  were  obtained  by  the  High-Speed  Flight 
Station  of  the  NACA  as  part  of  the  cooperative  program  with  the  Ames 
Aeronautical  Laboratory. 
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NOTATION 


A 

CL 

C1 

cn 

CY 

D 

% 

TZ 

iz 

Txz 

KZ 

KP/sa 

Kr/Sr 


aspect  ratio 

lift  coefficient 

rolling-moment  coefficient 

yawing-moment  coefficient 

side-force  coefficient 

differential  operator,  — 

dt 

moment  of  inertia  about  the  X  axis,  slug-ft2 

effective  rolling  moment  of  inertia  for  flexible  airplane, 
TX  -  q.QSbCz„,  slug-ft2 

moment  of  inertia  about  the  Z  axis ,  slug-ft2 

effective  yawing  moment  of  inertia  for  flexible  airplane, 
IZ  "  slug-ft2 

product  of  inertia,  slug-ft2 

dimensionless  radius  of  gyration  about  the  Z  axis,  wing 
spans 

gain  of  simplified  transfer  function 

°a 

gain  of  simplified  transfer  function 

Sr 


Kp/6r 

M 

S 

V 

w 

b 


gain  of  simplified  —  transfer  function 

5r 

Mach  number 
wing  area,  sq.  ft 
velocity,  ft/sec 
gross  weight  of  airplane,  lb 
wing  span,  ft 


MCA  TN  387^ 


local  chord  of  the  wing  measured  parallel  to  the  plane  of 
symmetry 

section  lift  coefficient 
acceleration  due  to  gravity,  ft/sec2 
mass  of  airplane,  slugs 

effective  mass  for  lateral  acceleration  of  flexible  airplane, 

normal  acceleration,  positive  downward,  gravity  units 

rolling  velocity,  radians /sec 

dynamic  pressure,  lb/sq  ft 

yawing  velocity,  radians /sec 

time,  sec 

structural  deflection,  positive  downward,  in. 
vertical  coordinate,  positive  downward,  in. 
angle  of  attack,  radians 

angle  of  sideslip,  radians  (except  as  otherwise  noted) 

total  aileron  deflection,  measured  in  a  plane  perpendicular 
to  the  hinge  line,  positive  right  aileron  up,  radians 
(except  as  otherwise  noted) 

rudder  deflection,  measured  in  a  plane  perpendicular  to  the 
hinge  line,  positive  trailing  edge  left,  radians  (except 
as  otherwise  noted) 

damping  ratio 

damping  ratio  of  numerator  term  of  —  transfer  function 

&a 

spanwise  coordinate,  in  wing  semispans 

parameter  used  in  evaluating  aerodynamic  derivatives  from 
transfer-function  coefficients,  defined  in  equation  (10) 

rolling-mode  time  constant,  sec 
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angle  of  bank,  radians 


§ 

/outputN 

\input/ 


phase  angle  of  output  quantity  minus  phase  angle  of  input 
quantity 


angle  of  yaw,  radians 
frequency,  radians /sec 

undamped  natural  frequency  of  numerator  term  of  ---  transfer 
function  a 


undamped  natural  frequency,  radians/sec 


C2p'  Cnp> 

S>  %’ 

V  V  °v 

xcp  ^  6cp 

c2..,  cn.. 

*  t 


derivative  of  coefficient  with  respect  to  subscript 


C C^,  C^,  derivative  of  coefficient  with  respect  to  subscript 

x  — 

2V 


-  CL  ,  per  sec 

2VIX^  V  ^ 


— r —  CL  ,  per  sec 
2VIY  T-t’ 

af 

q  Sb 

o  ^  2 

- -  CL  ,  per  sec 

\  V 


-  0^  ,  per  sec 

2 VI 7  np’ 


2Vl7_  “r 


Cn  ,  per  sec 


Nr 


— —  cn  ,  per  sec2 

17™  V 


q_  S 

_  C%,  per  see 


cp 


,  per  sec 


N, 


— —  C7  ,  per  sec2 

TXP  V 

40Sb  2 

- —  C  ,  per  sec 


Cj  ,  per  sec 


X, 


T  aSbC7 
IXZ  Ho  l 


-X- 


-xz 


Lp  +  rx  Np 


N. 


l0  +  %"» 

Nr  +  rZ-p^r 


V 

V 


Ne  +  Vp 

L6  +  rXFN8 


V 


N8  +  rZFL& 
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Subscripts 


B  body  axis 

eg  center  of 

t  tail 

wt  left  wing 

Dots  are  used  to 


example,  (3 


dp 

dt# 


gravity 

tip 

indicate  differentiation  with  respect  to  time;  for 


TEST  EQUIPMENT 


The  test  airplane  was  a  Boeing  B-k’JA  with  wing  vortex  generators 
installed  (fig.  l) .  An  external  nose  boom  and  an  optigraph  for  measuring 
the  movements  of  target  lights  mounted  on  the  wing  and  tail  were  added 
(figs.  1  and  2).  Pertinent  physical  characteristics  are  listed  in  table  I 
Left  and  right  aileron,  and  rudder  deflections  were  each  measured  by  three 
NACA  resistance-type  control-deflection  indicators  located  at  the  root, 
midspan,  and  tip.  An  examination  of  these  records  indicated  that  there 
was  a  negligible  difference  in  deflections  at  the  root,  midspan,  and  tip 
of  the  control  surfaces  for  the  flight  conditions  investigated.  Therefore 
only  the  midspan  deflection  records  were  used  in  the  analysis.  Wing  accel 
erations  were  measured  with  Statham  linear  accelerometers.  These  quanti¬ 
ties  were  recorded  on  Weston  12-channel  and  Consolidated  l8-channel  oscil¬ 
lographs.  Self-recording  turn  meters  were  used  to  measure  roll  and  yaw 
rates  at  the  center  of  gravity  and  tail.  Sideslip  angles  were  measured 
with  a  vane-type  pickup.  The  instruments  were  aligned  with  respect  to 
the  body  reference  line  and  locations  of  the  pertinent  instruments  are 
shown  in  figure  2. 


TEST  PROCEDURE 


Responses  to  both  aileron  and  rudder  pulses  were  obtained  for  a  Mach 
number  range  from  0.6  to  0.82  at  an  altitude  of  35,000  feet.  This  Mach 
number  range  is  that  in  which,  from  an  examination  of  wind-tunnel  data, 
the  aerodynamic  derivatives  could  be  expected  to  be  constant  for  a  given 
flight  condition,  and  thus  linear  analysis  techniques  could  be  used.  The 
center-of-gravity  location  was  approximately  21  percent  of  the  mean  aero¬ 
dynamic  chord  and  the  gross  weight  was  approximately  115,000  pounds.  The 
maneuvers  consisted  of  the  pilot  applying  a  sharp  aileron  or  rudder  pulse 
with  the  other  controls  held  fixed,  and  then  allowing  the  airplane  to 
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oscillate  for  about  25  seconds  (yaw  damper  off).  The  transient  motion 
did  not  completely  subside  in  this  time  because  of  the  lightly  damped 
Dutch  roll  mode.  Typical  time  histories  are  shown  in  figures  3  and 
In  order  to  show. the  structural  vibrations  more  clearly,  only  the  first 
portion  of  several  of  the  records  is  presented.  In  addition,  complete 
records  are  shown  for  roll  rate,  yaw  rate,  and  sideslip-angle  responses 
to  the  rudder  pulse  to  illustrate  the  motion  due  to  the  Dutch  roll  mode. 
For  some  of  the  maneuvers,  a  small  amount  of  Dutch  roll  motion  was  present 
when  the  control  input  was  applied.  This  effect  has  been  subtracted  out 
of  the  time  histories  shown  in  order  to  indicate  more  clearly  the  response 
to  the  control-surface  motions. 


PREDICTED  DYNAMIC  RESPONSES 


The  dynamic-response  characteristics  of  the  airplane  are  more  easily 
interpreted  if  the  time  histories  are  converted  to  frequency- response  or 
transfer-function  form.  In  order  to  establish  the  relationships  between 
these  forms  of  the  dynamic  response,  a  knowledge  of  the  equations  of 
motion  is  necessary.  These  equations  also  provide  a  basis  for  deter¬ 
mining  approximate  relations  for  transfer-function  coefficients  and 
aerodynamic  derivatives. 


Equations  of  Motion 


In  order  to  define  completely  the  motions  of  a  flexible  airplane,  a 
prohibitively  large  number  of  degrees  of  freedom  would  be  needed.  Hence, 
for  practical  purposes,  only  the  most  significant  degrees  of  freedom  can 
be  included.  A  simplification  of  particular  interest  from  the  standpoint 
of  control  studies  is  the  use  of  equations  of  motion  for  a  rigid  airplane 
with  the  coefficients  modified  to  include  the  pseudostatic  effects  of 
flexibility.  In  the  pseudostatic  method,  as  discussed  in  references  1 
and  2,  the  structural  deformations  are  assumed  to  be  in  phase  with  the 
aerodynamic  and  inertial  loads  resulting  from  airplane  motions.  The 
aerodynamic  forces  arising  from  these  deformations  are  included  by  modi¬ 
fying  the  derivatives  in  the  airplane  equations  of  motion. 

Three-degree- of -freedom  lateral- directional  equations  of  motion  with 
respect  to  stability  axes  are  given  in  Appendix  A.  For  the  flexible  air¬ 
plane  appropriate  derivatives  corresponding  to  structural  deflections  due 
to  inertial  and  dead-weight  forces  resulting  from  airplane  motions  have 
been  added  to  the  equations .  These  additional  terms  are  then  combined 
with  the  corresponding  rigid  airplane  inertias  to  form  an  effective  iner¬ 
tia  for  the  flexible  airplane.  As  explained  in  Appendix  A,  it  was  found 
to  be  permissible  to  neglect  a  number  of  the  inertial  derivatives  since  the 
cumulative  effect  of  these  quantities  on  the  transfer-function  coefficients 
was  found  to  be  negligible. 
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Aerodynamic  Derivatives 


The  stability  derivatives  were  obtained  by  available  theory  and  wind- 
tunnel  data.  Effects  of  structural  deflections  on  the  derivatives  were 
obtained  through  use  of  aerodynamic  and  structural  influence  coefficients 
in  the  manner  described  in  reference  1.  The  chordwise  centers  of  pressure 
of  all  aerodynamic  loads ,  except  loads  due  to  the  rudder  and  ailerons,  were 
assumed  to  be  at  the  quarter-chord  line  of  their  respective  surfaces.  The 
average  chordwise  center  of  pressure  for  loads  due  to  the  rudder  was  esti¬ 
mated  to  be  the  0A8  chord  from  the  pressure-distribution  data  given  in 
reference  5*  Although  no  pressure-distribution  data  were  available  for 
the  specific  wing-aileron  plan  form,  the  center  of  pressure  for  the  aileron 
loads  was  estimated  to  be  at  the  O.58  chord  from  an  examination  of  pressure 
distributions  on  similar  swept-wing  plan  forms.  Fuselage  influence  coeffi¬ 
cients  used  in  the  analysis  were  obtained  from  data  given  in  references  6 
and  7  and  are  listed  in  table  II  together  with  the  wing  influence  coeffi¬ 
cients  which  were  obtained  from  reference  8.  Wing  stations  and  equivalent 
weights  used  in  the  analysis  are  shown  in  figure  5*  An  example  of  the 
application  of  the  method  to  obtain  the  effect  of  wing  deflections  on  the 
derivative  C7  is  given  in  Appendix  B.  Predicted  aerodynamic  derivatives 

for  both  the  flexible  and  rigid  conditions  are  listed  in  table  III, 


Transfer-Function  Coefficients 


Transfer  functions  obtained  by  solving  the  three -degree -of- freedom 
equations  are  presented  in  Appendix  A.  As  noted  in  this  Appendix,  some 
of  the  transfer  functions  may  be  simplified  for  the  frequency  range  of 
interest.  After  the  spiral  mode  is  neglected,  rolling  response  to  aileron 
becomes 


p  __  A3  (D2  +  axD  4-  a2) 
ba  (D  +  Dr)(D2  +  ci D  +  c2) 


(1) 


The  transfer  function  may  also  be  written  as 


p_  S/sJ1  +  gga ?K  +  (pK)g] 

5a  (1  +  tD)[1  +  2pK  +  (D/u^)2] 


For  yawing  and  sideslipping  response  to  rudder  motion,  both  spiral  and 
rolling  modes  may  be  neglected  so  that 
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r  HpD 

5r  D2  +  CjD  +  c2 


(3) 


j.. _ 3» -  w 

5r  D2  +  c  jD  +  c2 


These  quantities  may  also  be  expressed  as 


r  V»rD 

8r  1  +  2^D/wn  +  (D/u^)2 


(5) 


J_  =  KP/5r _ 

5r  1  +  2p/un  +  (D/wn)2 


(6) 


The  predicted  frequency  responses  may  be  obtained  from  the  transfer 
functions  after  iu  is  substituted  for  the  operator  D.  The  conversion 
of  numerator  transfer-function  coefficients  from  stability  to  body  axes 
for  comparison  with  frequency  responses  measured  with  respect  to  body 
axes  is  given  in  Appendix  A. 


Approximate  Equations  to  Obtain  Aerodynamic  Derivatives 
From  Transfer-Function  Coefficients 


The  equations  for  the  aerodynamic  derivatives  were  obtained  as 
follows:  The  coefficients  C3,  C2,  and  Cx  of  the  characteristic  equa¬ 
tion  given  in  Appendix  A  were  expressed  in  terms  of  their  principal  aero¬ 
dynamic  derivatives.  The  coefficients  were  also  expressed  in  terms  of 
the  factors  of  the  characteristic  equation  Dr ,  c  x ,  and  c2.  By  equating 
the  coefficients  of  like  powers  of  D,  three  simultaneous  equations  were 
obtained  with  three  unknowns,  Lp,  Nr  +  Yp,  and  Npf.  A  cubic  equation 

for  Lp  was  then  obtained  by  eliminating  the  other  unknowns .  The  deriv¬ 
ative  Lp  was  found  by  using  an  approximate  solution  for  a  cubic  equation 
together  with  an  approximate  equation  for  the  amplitude  of  the  roll  to  yaw 
ratio  of  the  Dutch  roll  mode,  |£|.  Equations  for  the  two  other  derivatives 
could  then  be  obtained.  The  equations  are 
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S  =  _Dr  -  H 
Nr  +  "^p  =  |i  “  c1 

Np'  =  e2  +  M.(r>r  -  cx) 


where 


(7) 

(8) 
(9) 


(10) 


Note  that  these  equations  correspond  to  those  given  in  reference  k  except 
for  the  quantity  |i.  The  additional  term  (i  could  he  considered  as  a 
measure  of  the  coupling  between  the  Dutch  roll  and  rolling  modes.  The 
equation  used  for  Lp’  was  obtained  from  reference  k 


I  'J 1  + '  V7c2 


(11) 


Approximate  relations  for  control  effectiveness  were  obtained  by  substi¬ 
tuting  typical  numerical  values  in  the  equations  given  in  Appendix  A,  and 
the  following  expressions  (analogous  to  those  given  in  ref.  4)  were  found 
to  be  valid  for  the  present  airplane.  The  rudder  effectiveness  derivative 
was  obtained  from 


%r  =  Hr  or  -Jp 


(12) 


while  aileron  effectiveness  was  obtained  from 


Lsa  ■  A® 


(13) 


Aerodynamic  derivatives  may  be  found  from  the  quantities  on  the  left-hand 
side  of  equations  (7)  through  (13)  through  use  of  the  definitions  given  in 
the  notation. 


Thus  the  aerodynamic  derivatives  may  be  obtained  from  transfer- 
function  coefficients  matched  from  experimental  data,  provided  that  an 

estimate  is  made  of  the  product  of  inertia,  the  derivative  Cn  y  and  also 

3? 

the  effective  moments  of  inertia  for  the  flexible  airplane.  Estimated 
variations  in  effective  inertia  characteristics  are  presented  in  figure  6. 
The  principal  effect  of  flexibility  is  to  reduce  the  rolling  moment  of 
inertia.  This  effect  is  largely  due  to  the  mass  of  the  wing -mounted 
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nacelles.  In  order  to  demonstrate  the  adequacy  of  the  approximate 
equations ,  predicted  transfer-function  coefficients  were  first  calculated 
from  the  predicted  derivatives  (table  III)  using  the  complete  expressions 
for  the  transfer- function  coefficients.  Then  approximate  values  of  the 
derivatives  were  calculated  from  the  predicted  transfer-function  coeffi¬ 
cients  using  equations  (7)  to  (11),  and  also  using  the  equations  with 
|i  =  0,  which  correspond  to  the  equations  given  in  reference  4.  Comparisons 
of  these  derivatives  are  shown  in  figure  7.  Also  shown  are  derivatives 
calculated  from  equations  given  in  reference  4,  In  all  cases  the  values 
of  derivatives  using  the  approximate  equations  are  closer  to  the  estimated 
values  than  those  obtained  using  the  equations  from  reference  4.  This  is 
particularly  true  for  the  derivative ,  Cnr  +  ^K^Cy^.  Although  not  shown, 

a  similar  comparison  was  made  for  the  rigid  case,  and  the  inclusion  of  the 
quantity  \x  also  improved  the  agreement  with  the  predicted  derivatives. 


REDUCTION  OF  EXPERIMENTAL  DATA 
Frequency  Response 


The  Fourier  integral  was  used  to  convert  the  transient  responses  to 
frequency- response  form.  The  calculation  of  the  Fourier  integral  con¬ 
sisted  first  of  a  numerical  integration  over  the  first  part  of  the  record 
using  a  parabolic  fit  of  the  time  history  similar  to  that  described  in 
reference  9.  For  the  control  inputs  the  time  histories  showed  a  constant 
value  after  a  certain  time  and  hence  the  remaining  portion  of  the  integral 
was  evaluated  using  an  analytical  expression  for  the  small  step  displace¬ 
ment.  The  response  time  histories  did  not  reach  a  constant  value  since 
the  Dutch  roll  oscillation  was  very  lightly  damped,  and  analytical  expres¬ 
sions  for  both  a  step  and  a  damped  sinusoidal  oscillation  (e.g. ,  ref.  10) 
were  used  to  evaluate  the  remaining  part  of  the  integrals.  These  calcu¬ 
lations  were  performed  on  a  card- programmed  IBM  650  digital  computer. 

All  traces  were  read  at  0.05-second  intervals  except  sideslip  angle, 
which  was  read  at  0.1-second  intervals.  In  some  cases,  Fourier  transforms 
were  obtained  from  roll  and  yaw  angular-acceleration  traces  in  order  to 
define  better  the  corresponding  angular  rate  traces  at  higher  frequencies. 
Results  from  sideslip  records  at  the  higher  frequencies  are  not  presented 
because  of  the  lack  of  well-defined  trends.  This  is  probably  due  to  the 
relative  attenuation  of  the  higher  frequency  amplitudes  measured  by  an 
angular  displacement  instrument  as  compared  with  those  measured  by  angular 
rate  instruments.  Phase-angle  corrections  were  made  for  the  dynamic 
response  of  the  instruments. 
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Transfer-Function  Coefficients  and 
Aerodynamic  Derivatives 


Transfer-function  coefficients  of  the  simplified  transfer  functions 
for  p/5a,  r/&r,  an<3*  3/^r  described  in  a  previous  section  were  obtained 
by  matching  the  response  time  histories  (e.g.,  figs.  3  and  4)  by  means  of 
an  analog  computer.  They  were  then  converted  from  body  axes  to  stability 
axes.  In  some  cases ,  small  residual  oscillations,  which  were  principally 
due  to  the  lightly  damped  Dutch  roll  mode,  were  present  when  the  airplane 
control  pulse  was  applied.  Initial  conditions  were  applied  to  the  simu¬ 
lator  to  include  this  motion,  although  in  all  cases  this  had  only  a  small 
effect  on  transfer-function  coefficients.  In  matching  the  time  histories 
of  p/ba  on  the  analog  computer,  it  was  found  convenient  to  represent  the 
third-order  transfer  function  by  parallel  networks  of  first-order  and 
second-order  transfer  functions.  This  representation  is  equivalent  to 
separating  the  third-order  transfer  function  into  the  partial  fraction 
form  shown  in  Appendix  A.  In  this  way  contributions  of  the  Dutch  roll 
and  rolling  modes  to  the  resultant  motions  could  be  matched  separately. 

Aerodynamic  derivatives  were  then  obtained  from  the  appropriate 
transfer- function  coefficients  using  the  approximate  equations  (eqs.  (7) 
to  (13))  given  in  the  previous  section. 


RESULTS  AM)  DISCUSSION 
Frequency  Responses 


In  the  discussion  that  follows,  effects  of  the  relatively  low 
frequency  lateral- directional  modes  will  be  investigated  first.  Com¬ 
parisons  will  be  made  between  measured  and  predicted  responses  to  deter¬ 
mine  the  frequency  range  for  which  the  predictions  are  adequate.  Then 
the  higher  frequency  structural  modes  will  be  discussed  and  compared  with 
values  obtained  from  the  ground  oscillation  tests  of  reference  11.  The 
frequency- response  data  to  be  presented  include:  measured  frequency 
responses  to  aileron  and  rudder  inputs  over  a  fairly  wide  Mach  number 
range  (figs.  8  and  9) >  typical  comparisons  of  measured  responses  at  the 
center  of  gravity  and  tail  (figs.  10  and  11),  and  typical  comparisons  of 
measured  and  estimated  frequency  responses  (figs.  12  and  13). 

Lateral-directional  modes . -  The  ailerons  excite  principally  the 
rolling  mode  as  well  as  some  Dutch  roll  mode  as  shown  in  the  roll- rate 
response  (fig.  8(a)),  which  is  seen  to  be  of  the  same  form  as  that  indi¬ 
cated  in  equation  (2).  The  rolling  mode  (a  first-order  term)  causes  the 
net  shift  of  90°  in  the  phase  angle  between  the  frequencies  of  0.5  to  5 
radians  per  second  and  the  slope  of  the  amplitude  curve  in  the  frequency 
range  of  2  to  7  radians  per  second,  while  the  peak  in  the  amplitude  curve 
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that  occurs  at  frequencies  from  about  1  to  1.5  radians  per  second  is  due 
to  the  Dutch  roll  mode  (a  second-order  term) .  Responses  in  sideslip  and 
yaw  rate  are  quite  small  (figs.  8(b)  and  8(c)).  The  rudder  mainly  excites 
the  Dutch  roll  mode  (fig.  9).  For  both  inputs,  differences  in  tail  and 
center- of -gravity  roll  and  yaw  rates  are  small  in  this  frequency  range. 
This  may  also  be  seen  from  the  time  histories  (figs.  3  and  *0  • 


The  predicted  results  for  the  flexible  airplane  in  general  agree 
reasonably  well  with  experimental  results  for  frequencies  up  to  about  3 
to  8  radians  per  second,  with  the  possible  exception  of  wing-tip  accel¬ 
eration  (figs.  12(e)  and  13(f)).  Also  some  discrepancies  are  apparent 
in  the  lowest  frequency  range  (near  0.5  radian/sec)  which  is  the  region 
in  which  the  experimental  frequency-response  accuracy  obtained  from  pulse' 
type  inputs  becomes  relatively  poor,  as  was  discussed  in  references  1 
and  4.  For  the  aileron  responses,  the  predicted  gain  for  the  flexible 
case  is  much  lower  than  that  for  the  rigid,  reflecting  the  reduction  in 
aileron  effectiveness.  However,  the  p/5r  responses  (fig.  13(a))  indi¬ 
cate  a  higher  gain  for  the  flexible  case.  This  is  largely  due  to  the 
higher  value  of  the  ratio  of  to  the  effective  moment  of  inertia  in 


the  flexible  case.  Flexibility  also  reduces  the  frequency  and  damping  of 
the  Dutch  roll  mode.  A  discrepancy  between  theory  and  experiment  in  the 
curves  for  r/6a  (fig.  12(b))  is  that  the  experimental  phase  angles 
shift  l80°  while  the  predicted  values  approach  a  3^0°  shift.  One  possi¬ 
ble  explanation  is  that  the  smaller  phase  shift  could  have  been  obtained 
from  the  predicted  curve  if  a  more  negative  value  of  Cn  had  been 

^a 

assumed.  A  structural  mode  could  also  affect  the  phase  angle  in  this 
region. 


Structural  modes.-  A  peak  occurs  in  all  the  frequency  responses  to 
aileron  (fig.  8)  at  about  13.4  radians  per  second  which  will  be  called 
the  first  wing  antisymmetric  bending  mode.  This  compares  with  l4.4 
radians  per  second  obtained  from  ground  vibration  tests  (ref.  11).  The 
mode  is  particularly  noticeable  in  the  center- of -gravity  roll  rate  and 
wing- tip  acceleration  responses.  In  the  time  history  for  the  wing  opti- 
graph,  this  mode  is  seen  to  predominate.  Another  point  of  interest  is 
that  the  amplitude  of  the  yaw  rate  at  the  tail  is  greater  than  that  for 
the  center  of  gravity,  as  may  be  seen  from  both  the  frequency- response 
and  time-history  curves.  The  mode  frequency  remains  essentially  constant 
for  the  flight  range  investigated  although  there  is  a  slight  increase  in 
amplitude  with  Mach  number.  Measurements  of  the  time  history  indicated 
that  the  damping  was  also  essentially  constant  over  the  flight  range  at 
a  value  of  about  0.025  critical. 


The  shape  of  this  mode,  obtained  from  measurements  of  the  free 
oscillation  amplitudes  of  the  accelerometer  and  turn-meter  traces  result¬ 
ing  from  aileron  inputs,  is  shown  in  figure  l4.  When  the  curve  was  deter¬ 
mined,  small  corrections  for  chordwise  displacement  of  the  instruments  were 
made  where  necessary  to  transfer  the  accelerometer  results  to  the  17- 
percent  and  58-percent  chord  lines  corresponding  to  the  front  and  rear 


NACA  TN  387^ 


15 


spars,  respectively.  In  the  placement  of  pickups  for  an  autopilot  system, 
it  is  desirable  to  locate  them  at  points  on  the  structure  that  are  least 
affected  by  the  structural  modes  in  order  to  avoid  introducing  undesired 
signals.  For  a  roll  gyro,  the  optimum  location  for  minimizing  the  angular 
velocity  input  of  the  antisymmetric  wing  mode  would  be  where  the  wing  mode 
shape  has  zero  slope  which,  as  shown  in  figure  l4,  is  at  rj  =  O.38  near 
the  inboard  nacelle,  while  for  a  linear  accelerometer,  the  optimum  loca¬ 
tion  would  be  at  a  nodal  point,  which  in  figure  l4  is  at  T)  of  about  0.55 
just  outboard  of  the  nacelle.  The  wing  accelerometer  node  line  obtained 
from  the  two  span  positions  shows  good  agreement  with  that  obtained  from 
ground  vibration  tests. 

In  the  wing- tip  acceleration  curve  (fig.  8(d)),  a  dip  in  the  amplitude 
and  a  shift  of  approximately  l80°  in  the  phase  angle  relative  to  the 
center-of-gravity  roll  rate  (fig.  8(a))  occurs  in  a  frequency  range  some¬ 
what  lower  than  the  wing  bending  frequency.  This  motion  can  be  illustrated 
as  follows:  At  low  frequencies,  the  relation  between  wing  bending  and  roll 
angle  is  as  shown  in  sketch  (a).  As  the  frequency  increases,  the  amount  of 


wing  bending  relative  to  roll  angle  becomes  greater  so  that  the  accelera¬ 
tion  at  the  tip  reverses  sign  with  respect  to  the  roll  angle  (sketch  (b)). 


This  trend  continues  until  the  bending  frequency  is  reached. 

An  examination  of  figure  8  shows  that  no  additional  modes  were  excited 
to  any  appreciable  extent  in  the  measured  frequency  range  above  that  of  the 
wing  first  antisymmetric  bending  mode,  although  two  modes  were  measured 
from  ground  vibration  tests  (ref.  11). 

While  the  rudder  pulses  also  excite  the  wing  antisymmetric  bending 
mode,  their  principal  effect  is  to  excite  two  slightly  higher  frequency 
modes  which  are  predominantly  fuselage  modes  as  may  be  seen  from  the 
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responses  of  yaw  and  roll  rate  at  the  tail  and  center  of  gravity  (figs.  9 
and  11) .  The  lower  frequency,  which  is  predominantly  fuselage  side  bend¬ 
ing,  is  16.5  radians  per  second  while  the  upper  frequency,  which  is  prin¬ 
cipally  fuselage  torsion,  is  23  radians  per  second.  These  frequencies 
compare  with  ground  vibration  test  values  of  18.0  and  24.5  radians  per 
second,  respectively  (ref.  11).  The  prominance  of  these  modes  at  the 
tail  may  also  be  seen  from  the  transient  responses  due  to  the  rudder 
input  (fig.  4). 

The  frequency  responses  for  roll  and  yaw  rate  at  the  tail  due  to  a 
rudder  pulse  (fig.  11)  also  indicate  regions  of  low  amplitude  at  frequen¬ 
cies  below  the  structural  natural  frequencies  (about  3  to  8  radians  per 
second) .  While  phase  angle  for  the  tail  yaw  rate  shifts  l80°  out  of  phase 
with  that  of  the  center  of  gravity,  as  was  the  case  for  the  wing  at  fre¬ 
quencies  below  that  of  the  antisymmetric  bending  mode,  the  phase  angle 
shifts  for  the  tail  and  center- of -gravity  roll  rate  are  approximately  the 
same  so  that  they  remain  in  phase.  It  is  of  interest  to  note  that  by  the 
inclusion  of  only  the  pseudostatic  effects  of  sideslip,  yawing  accelera¬ 
tion,  etc.,  on  the  tail  deflections,  these  amplitude  and  phase-angle  trends 
are  predicted  (figs.  13(b)  and  13(d)).  While  the  direction  of  the  l80° 
phase-angle  shift  does  not  agree  for  the  tail  yaw  rate,  the  resultant  shift 
is  the  same.  Since  the  amplitude  ratio  is  quite  low  in  this  region,  some 
question  does  exist  as  to  the  direction  of  the  phase-angle  shifts  for  both 
the  experimental  and  predicted  values. 


Transfer-Function  Coefficients 


Comparisons  of  experimentally  determined  transfer-function  coeffi¬ 
cients  with  predicted  values  (with  respect  to  stability  axes)  for  both 
the  flexible  and  rigid  case  based  on  the  simplified  transfer  functions 
(eqs.  (2),  (5),  and  (6))  are  shown  in  figure  15.  In  general,  agreement 
between  experiment  and  theory  for  both  numerator  and  denominator  is  seen 
to  be  good,  as  was  the  case  reported  in  reference  2  for  25,000  feet. 

For  the  denominator  terms ,  the  values  of  measured  damping  ratio  are 
an  average  of  0.015  less  than  the  estimated  flexible  values.  While  the 
ratios  of  measured  to  estimated  values  are  relatively  large,  the  average 
magnitude  is  approaching  zero  at  the  lower  values  of  Mach  number,  and 
hence  the  discrepancy  is  considered  to  be  small.  Also  some  scatter 
occurred  in  the  time  constant  t.  This  parameter  was  difficult  to  eval¬ 
uate  because  it  did  not  have  as  important  an  effect  on  the  transient 
response  to  a  pulse  as  did  the  other  parameters.  Flexibility  is  seen  to 
reduce  both  the  frequency  and  damping  of  the  Dutch  roll  mode.  These 
reductions  are  caused  primarily  by  the  decrease  in  vertical-tail  effec¬ 
tiveness  due  to  fuselage  and  tail  flexibility  having  a  larger  effect  than 
the  decrease  in  effective  yawing  moment  of  inertia.  Flexibility  changes 
the  time  constant  t  very  little.  This  term  depends  principally  on  the 
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ratio  of  effective  rolling  moment  of  inertia  to  the  aerodynamic  damping 
in  roll.  Since  the  rolling  moment  of  inertia  is  due  largely  to  the 
nacelle  and  wing  weights,  wing  flexibility  reduces  both  the  aerodynamic 
and  inertia  loads  by  about  the  same  amount. 

For  the  numerator  terms  also,  good  agreement  is  obtained  between 
theory  and  experiment.  The  scatter  that  occurred  in  the  measured  values 
of  £a  resulted  from  the  small  effect  of  £a  on  the  time  histories 
matched  with  the  analog  computer.  The  comparison  of  the  matched  curves 
in  figures  3  and  k  with  the  measured  time  histories  indicates  that  the 
transfer  functions  corresponding  to  the  Dutch  roll  mode  only  for  the  rud¬ 
der  pulses,  and  the  rolling  and  Dutch  roll  modes  for  the  aileron  pulses 
were  sufficient  to  define  the  time  histories  adequately.  The  greatest 
effect  of  flexibility  is  seen  to  be  in  the  reduction  in  gain  for  p/6a. 
This  reflects  the  reduction  in  aileron  effectiveness  associated  with  out¬ 
board  ailerons.  The  predicted  values  of  £a  and  wa  are  slightly  differ¬ 
ent  from  £  and  u)n  and  approach  these  quantities  at  the  higher  Mach  num¬ 
bers.  These  differences  depend  upon  the  aerodynamic  derivatives  that 
effect  the  quantities  q  and  Cn  ,  which  become  relatively  small  at  the 

^a 

higher  Mach  numbers  and  corresponding  low  angles  of  attack.  These  effects 
are  reflected  in  the  time  histories  by  the  amount  of  excitation  of  the 
Dutch  roll  mode  and  also  by  a  reduction  in  over-all  gain.  For  the  limit¬ 
ing  case  of  £a  =  £  and  wa  =  effects  of  the  Dutch  roll  mode  are  elim¬ 
inated  and  the  transfer  function  (eq.  (2))  reduces  to  a  first-order  form. 


Aerodynamic  Derivatives 


Generally  good  agreement  is  obtained  between  predicted  and  experi¬ 
mental  values  of  the  aerodynamic  derivatives  (fig.  l6) .  Differences 
between  predicted  aerodynamic  derivatives  for  the  flexible  and  rigid 
cases  are  generally  somewhat  greater  than  the  differences  between 
transfer-function  coefficients  since  the  latter  are  partially  compensated 
by  changes  in  the  effective  moments  of  inertia  due  to  flexibility.  In 
general,  flexibility  tends  to  reduce  the  values  of  the  aerodynamic  deriv¬ 
atives  with  the  largest  reduction  occurring  for  C^  and  C^  .  However, 

^a  2 

flexibility  tends  to  increase  the  magnitude  of  C7  since  the  usual 

decrease  in  value  of  a  derivative  due  to  flexibility  is  more  than  compen¬ 
sated  by  the  increase  in  dihedral  due  to  wing  lift. 


CONCLUSIONS 


The  dynamic  lateral- directional  responses  to  aileron  and  rudder  pulses 
of  a  large  flexible  airplane  for  flight  conditions  of  0.6  to  0.82  Mach  num¬ 
ber  and  an  altitude  of  35^000  feet  have  been  investigated  and  compared  with 
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predicted  responses  leased  on  the  lateral-directional  equations  of  motion. 
Pseudostatic  effects  of  flexibility  were  included  through  use  of  struc¬ 
tural  and  aerodynamic  influence  coefficients.  The  following  conclusions 
have  been  made: 

1.  The  manually  applied  pulse  technique  was  found  to  be  adequate  to 
obtain  data  suitable  for  defining  frequency  responses  for  the  frequency 
range  from  1  to  25  radians  per  second.  This  range  includes  the  relatively 
low  frequency  Dutch  roll  and  rolling  modes,  as  well  as  the  higher  frequency 
wing  first  antisymmetric  bending,  fuselage  side-bending,  and  fuselage 
torsion  modes. 

2.  Reasonably  good  agreement  was  obtained  in  the  frequency  range 
below  about  one  half  of  the  lowest  structural  frequency  between  measured 
and  estimated  frequency  responses. 

3.  The  principal  effect  of  flexibility  on  predicted  transfer-function 
coefficients  is  the  reduction  in  gain  of  airplane  response  to  aileron 
motion,  reflecting  the  reduction  in  aileron  effectiveness,  and  also  a 
reduction  in  damping  and  frequency  of  the  Dutch  roll  mode. 

4.  The  principal  lateral-directional  derivatives  and  control- 
effectiveness  derivatives,  which  were  evaluated  from  measured  transfer 
functions  through  use  of  approximate  equations,  indicated  good  agreement, 
in  general,  with  predicted  derivatives  when  the  pseudostatic  effects  of 
flexibility  were  included. 

5.  In  the  frequency  range  above  half  the  lowest  structural  frequency, 
the  effects  of  structural  modes  were  sufficiently  large  (even  at  the  cen¬ 
ter  of  gravity)  for  the  responses  measured,  so  that  predictions  of  the  fre¬ 
quency  responses,  which  included  only  the  pseudostatic  effects  of  flexi¬ 
bility,  were  inadequate. 

6.  The  frequencies  of  the  three  structural  modes  measured  in  flight 
were  within  about  10  percent  of  the  frequencies  of  the  corresponding  modes 
measured  in  ground  vibration  tests.  In  addition,  the  wing  node  line  for 
the  wing  first  antisymmetric  bending  mode  measured  in  flight  showed  good 
agreement  with  that  obtained  from  the  ground  vibration  test. 
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APPENDIX  A 

PREDICTED  AIRPLANE  RESPONSE 
Equations  of  Motion 


The  three  lateral-directional  equations  of  motion  with  respect  to 
stability  axes  given  in  Appendix  C  of  reference  4  for  a  rigid  airplane 
will  be  modified  for  use  for  a  flexible  airplane.  For  the  flexible  case, 
additional  terms  must  be  added  to  take  account  of  aerodynamic  forces 

resulting  from  structural  deflections  due  to  inertial  loads  (C-,.. ,  C-?.., 

^cp  L\ 

Cn..,  Cy#,  and  Cy#)  and  dead-weight  loads  (Cy  ) .  A  more  complete  evalua- 

♦  ♦  0  T 

tion  of  these  additional  derivatives  would  also  have  included  other  deriv¬ 
atives  dependent  on  cp,  \|?,  \Jr+(3,  and  cp.  However,  by  the  substitution  of 
typical  numerical  values  in  the  transfer-function  coefficients,  the  cumu¬ 
lative  effects  of  these  other  derivatives  on  the  transfer-function  coef¬ 
ficients  were  found  to  be  negligible  even  though  some  of  them  were  of 

appreciable  size  (e.g.,  Cn.  and  Cn.  were  several  times  larger  than 

\|/  0 

Cn-rCV^V)  for  the  largest  values  of  M) .  The  rolling  moment,  yawing 

moment,  and  side-force  equations,  with  flight-path  angle  assumed  zero, 
may  be  expressed  as 


IVD2  -  <1  SbC7  D2  -  q.  SbC7 
a  o  o  i 


P  2V 


D  )  cp  + 


-JXZP  ~  VbCZvD  ‘  '1oSbCZ.  iv  * 


*oSby  =  Vbcz&8  <A1) 


(■ 


-IXZd2  “  *oSbCnp  2y  D  ^  Cf)  + 


IZD2  -  q  SbCn>D2  -  q  SbC  ^  D  U  -  a  SbCn  B  =  a  SbCn  8 
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The  remaining  inertial  deflection  derivatives  were  then  combined  with  the 
mass  terms  to  form  effective  inertias  for  the  flexible  airplane.  Thus 


%  '  h  '  qosi% 


mF  =  m  - 


q  S 

—  Cy  ♦ 
v  Yp 


<loS  . 

or  m - Cv 

V  1 T 

Since  the  derivatives  CY  and  CY  each  represent  aerodynamic  forces  due 

i  p 

to  lateral  acceleration,  they  are  of  equal  magnitude,  and  hence  rap  may 
he  expressed  in  terms  of  either  derivative.  The  three  equations  (Al) , 
(A2),  and  (A3)  can  then  he  written  in  a  more  convenient  form  hy  dividing 
them  hy  Ix  ,  I^,  and  m^V,  respectively,  and  introducing  new  symbols 


(d2  -  LpD^  cp  +  (-r^D2  -  i  “  L^p  =  L&S 


(A4) 


(-r^D2  -  Npd)  cp  +  (d2  -  i  -  Npp  =  Ng8 


(A5) 


-Y^cp  +  D+  +  (D  -  Yp)p  =  Y6S 


(A 6) 


Note  that  the  term  (which  equals  g/v)  remains  the  same  as  for  the 
rigid  case  since  the  derivative  Cy^  is  due  to  dead-weight  loads  that 
are  distributed  in  the  same  manner  as  the  lateral  acceleration  loads. 


Transfer  Functions 


Transfer- function  coefficients  for  cp/& ^  and  3/6  can  be  obtained 

by  simultaneous  solution  of  the  three  equations  (kb) ,  (A5) ,  and  (a6) . 
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Denominator  coefficients The  denominator  of  the  transfer  functions, 
A,  may  be  expressed  as 

A  =  D(C4D4  +  C3D3  +  C2D2  +  CjD  +  CQ) 

where 


C2  =  V  +  Yp(Lp'  +  Nr»)  +  LpNr  -  WpLj. 

Ci  =  LPNP  '  V*  '  W  +  Yp(V*p  -  NrLp) 

Co  =  Ycp(LpNr  -  NpLp) 


For  the  moderate  angle- of- attack  range  considered,  rv  and 


are 


small  quantities  and  hence 
in  factored  form  as 


C4  ~  1 


The  denominator  can  then  be  expressed 


A  =  D(D  +  DS)(D  +  Dr)(D2  +  cxD  +  c2) 

where  Ds  and  Dr  represent  the  spiral  and  rolling  modes,  respectively, 
and  cx  and  c2  are  coefficients  that  define  the  damping  and  frequency 
of  the  Dutch  roll  mode. 


Numerator  coefficients .  -  With  p  =  Dcp,  and  r  =  D^,  the  numerator 
coefficients  of  the  following  transfer  functions  can  be  expressed  in 
terms  of  the  derivatives. 


p  _  DtAgD3  +  AjjD2  +  AXD) 
8  "  A 


where 


A2  =  -YpL6,  +  N5Lr  "  LBNr  +  Y6Lp' 

Ai  “  ^6^3  +  “  N5^r)  +  Yg(^r^p  ” 
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r  _  D(B3D3  +  BgD2  +  BXD  +  B0) 

"b  a 

where 

B3  =  V 

B2  =  L8NP  "  NSLP  ”  Y(3N5  +  Y8Np 

Bi  =  W*  ”  L&NP^  +  Y5^LpNP  “  NPLP^ 

Bo  =  Wp  “  N&V 

p  _  D(EsD3  +  EpD2  +  EjD  +  Ep) 

5  ~  A 


where 

E3  =  "^bC1  "  rxFrzF) 

E2  =  -Yg(Nr*  +  Lp*)  “  Ng* 

Ex  =  YqjLg*  +  NgLp  -  LgNp  +  Yg(lpNr  -  NpLj.) 

Eo  =  W*  -  LSNr) 


Simplifications  of  the  transfer  functions.-  When  transfer  functions 
are  evaluated  from  measured  data,  it  is  desirable  to  use  as  simple  a  form 
as  possible  which  will  still  adequately  fit  the  data.  Calculations  indi¬ 
cated  that  the  spiral  mode  factor  Ds  was  very  small  and  could  be  neg¬ 
lected  for  the  frequency  range  of  interest.  Thus  the  rolling  response  to 
aileron  was  simplified  to 

p  _  A3(D2  +  axD  +  ag) 

8a  (D  +  Dr)(D2  +  CjD  +  c2) 


where 
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To  determine  estimated  values  for  use  with  the  curve  fitting  of  the 
measured  responses  on  the  analog  computer ,  a  partial  fraction  form  of 
the  transfer  function  was  advantageous  to  use. 


D  +  Bp 


HpD  +  J 
B2  +  CpD  +  c2 


In  obtaining  approximate  relations  for  r/br  and  p/6r,  it  was 
desirable  first  to  write  the  transfer  functions  in  partial  fraction 
form. 


'p  Gp  HpB  +  Jp 

+  -  + 


5r  B  +  Bs  B  +  Bp  B2  +  CpD  +  c. 


p  fp 


Gp  HpD + Jp 


B  +  Bq  B  +  B^»  B2  +  c  2.B  +  C2 


By  the  substitution  of  typical  numerical  values,  all  numerator  terms  were 
found  to  be  negligible  except  Hr  and  J^.  Thus ,  the  transfer- function 
coefficients  could  be  simplified  to 


p  HpB 

&r  D2  +  cxD  +  c2 

_P_  =  JP 

&r  D2  +  cxD  +  c2 


2k 
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Transfer  functions  for  the  rear  part  of  the  fuselage  and  the  wing 
tip.-  In  addition  to  the  transfer-function  coefficients  for  quantities 
at  the  center  of  gravity,  the  coefficients  for  r^/8,  P^/h >  and  n^/8 

needed.  The  evaluation  of  these  quantities  includes  responses  at  the 
center  of  gravity  and  also  responses  due  to  pseudostatic  structural 
deflections  at  the  particular  location.  The  equation  for  yaw  rate  at 
the  rear  part  of  the  fuselage  is 


are 


r  DA\|^ 
5  +  8 


where  A\|r  represents  the  total  change  in  angle  of  yaw  in  radians  at  the 
t 

rear  part  of  the  fuselage  due  to  structural  deformations  resulting  from 
a  control  input.  The  following  quantities  were  included  in  determining 


^t  A^t  p  t  r  |  Att  j  |  A^t 

8  (38  r8  $  8  8 


at  the  rear  part  of  the  fuselage  due  to  structural  deformations  resulting 
from  a  unit  change  in  p. 


Similarly,  the  equation  for  roll  rate  at  the  rear  part  of  the  fuselage 
is 

ft  _  p  f^t 

5  5 +  5 

where  Acp  represents  the  total  change  in  angle  of  roll  at  the  rear  part 
t 

of  the  fuselage  due  to  structural  deformations  resulting  from  a  control 
input  and  may  he  expressed  as 

5  p  8  r  &  t  8  5 


MCA  TN  387^ 


25 


The  equation  for  acceleration  at  the  left  wing  tip  is 

=  _  Dp  D2Zvt 

5  8  125 

Calculations  indicated  that  the  effect  of  D2z  ,/b  was  small  for  the 

wt' 

frequency  range  up  to  the  first  antisymmetric  bending  frequency  and, 
hence,  was  neglected. 

Conversion  of  transfer-f unction  coefficients  from  stability  axes  to 
body  axes.-  In  order  to  compare  predicted  results  with  measured  frequency 
responses  which  were  obtained  with  respect  to  body  axes,  the  predicted 
transfer  functions  were  converted  from  stability-axis  to  body-axis  form. 
Since,  for  the  angle- of- attack  range  considered,  the  quantity  C4  was 
essentially  unity,  only  the  numerator  terms  of  the  transfer  function  need 
be  modified.  The  equations  for  conversion  from  stability  axes  to  body 
axes  are  (e.g.,  ref.  4) 


PB 


=  p  cos  a  -  r  sin  a 


r^  -  r  cos  a  +  p  sin  a 


tan  PB 


tan  3 
cos  a 


For  small  angles  of  attack,  the  equations  may  be  simplified  to 


PB  =  p  -  ra 
rB  =  r  +  pa 

eB  =  p 


These  relations  will  also  hold  for  derivatives  of  the  angles.  Hence,  the 
numerator  coefficients  for  p/5  and  r/6  may  be  converted  as  follows: 


Bfa 


—  B^  +  Aj_a 
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APPENDIX  B 


EXAMPLE  CALCULATION  OF  WING  CONTRIBUTION  TO  C7  USING 
AERODYNAMIC  AND  STRUCTURAL  INFLUENCE  COEFFICIENTS 


A  general  description  of  the  method  used  as  applied  to  symmetrical 
loadings  was  given  in  reference  1.  An  example  application  to  obtain 
antisymmetric  loadings  will  he  presented  here.  The  equations  to  be  used 
will  be  given  first,  and  then  a  numerical  example  to  obtain  C 1  will 
be  shown.  p 

The  rolling-moment  coefficient  resulting  from  an  arbitrary  antisym¬ 
metric  angle -of- attack  distribution  may  be  expressed  as  a  power  series 
of  the  dynamic  pressure  q. 


Cl 


F 


+  + 

lA 


V2  +  v3  + 


(Bl) 


where 

rolling-moment  coefficient  for  the  flexible  wing 


C7  rolling-moment  coefficient  for  the  rigid  wing 

*R 


Cj  q  increment  in  rolling-moment  coefficient  resulting  from  structural 
A  deflections  due  to  the  Cj  loading 


ci^2 

LB 


increment  in  rolling-moment  coefficient  resulting  from  structural 

deflections  due  to  the  C^  loading 

A 


For  the  wing  structure  considered  in  the  present  example,  the  ratios 
of  all  terms  after  the  first  one.  -C7  /C7  ,  -C7  /C7  ,  etc.,  are  essen- 

lA  lB  lB  lC 

tially  equal  to  a  single  constant  k.  This  result  may  be  interpreted  to 
mean  that  the  successive  wing-deflection  shapes  that  produce  these  loads 
are  essentially  the  same,  and,  hence,  this  portion  of  the  wing  loading 
may  be  represented  as  a  single-degree-of-freedom  system.  Thus  the  incre¬ 
mental  rolling-moment  coefficient  due  to  this  portion  of  the  wing  deflec¬ 
tions  may  be  represented  by  the  expression  Cj  q/l+kq  and  it  will  be 

valid  for  large  positive  values  of  the  quantity  kq  even  though  the  series 
would  not  converge  for  positive  kq  greater  than  1.  The  equation  for 
rolling-moment  coefficient  for  the  flexible  wing  will  then  be 
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C 


l 


F 


1  +  kq 


(B2) 


The  quantities  in  equation  (B2)  will  now  be  determined  through  use  of 
aerodynamic  and  structural  influence  coefficients . 

Antisymmetric  aerodynamic  influence  coefficients  were  obtained  through 

use  of  reference  12.  First,  p  angle -of- -attack  influence  coefficients, 

7  vn  0  7 

which  represent  the  angle  of  attack  at  station  v  due  to  a  continuous 
loading  function  having  a  unit  value  at  station  n  and  zero  value  at  the 
other  stations,  were  obtained  from  reference  12.  These  were  then  con¬ 
verted  to  loading  influence  coefficients  Avn  (i.e.,  the  loading  c  c/2b, 

at  station,  n,  due  to  a  continuous  angle-of-attack  distribution  with  a 
unit  value  at  station  v  and  zero  values  at  the  other  stations)  by  solving 
for  the  span  loadings  for  a  unit  angle  of  attack  at  one  control  station. 
This  was  done  for  each  of  the  spanwise  control  stations  1,  2,  and  3  located 
at  T)  =  0.92^,  0.707,  and  O.383,  respectively.  The  resulting  rolling- 
moment  coefficient  for  a  unit  angle  of  attack  at  a  control  station,  v,  can 
then  be  calculated  using  the  following  equation,  which  has  the  form  of 
equation  15  in  reference  12. 


C 


ly 


=  M  [A  , 

16  1  V2 


0.707(AV1  +  AV3)] 


(B3) 


The  wing  structural  influence  coefficients  are  listed  in  table  II. 
For  use  with  the  aerodynamic  influence  coefficients,  they  were  converted 
to  the  form  of  a  change  in  streamwise  angle  of  attack  in  radians,  due  to 
a  1000-pound  load  at  a  given  percent  chord  position.  All  aerodynamic 
loads  were  assumed  to  act  along  the  quarter- chord  line  except  for  loads 
due  to  the  ailerons.  The  change  in  streamwise  angle  of  attack  at  sta¬ 
tion  m  due  to  the  loading  resulting  from  an  angle  of  attack  at  station 
v,  with  q  =  1  psi,  may  be  expressed  as 

(B4) 


(O 


vm 


=  b2(!2) 
1000 


I 


n=i 


^vnBmn^n 


with 


m  =  1,  2,  3 

v  =  1,  2,  3 
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where 

Avn  aerodynamic  influence  coefficients  for  the  loading  function  at 
station  n  due  to  an  angle  of  attack  at  station  v 

Smn  structural  influence  coefficients  expressing  the  change  in  angle 
of  attack  at  station  m  in  radians  due  to  a  1000-pound  load 
at  the  quarter- chord  station  n 

In  integrating  factors  with  the  following  values 

I-L  =  0.1502 

12  =  0.2776 

13  =  0.3628 

The  aerodynamic  loading  due  to  (en)  is  then  obtained  from  the  summation 

u  vm 

3 

(A/0  vn  =  ^  (eo)vnAnn  ^5) 

m=i 


from  which  the  incremental  rolling-moment  coefficient  may  be  obtained. 


(%)  ■  ft  {w.  +  °-707[(aa>v1  +  uoj  } 


(b6) 


The  next  incremental  rolling-moment  coefficient 

resulting  from  the  loading  coefficients  (A&) vn 
manner . 


b)v 


due  to  deflection 


is  calculated  in  a  similar 


Thus,  the  rolling-moment  coefficients  resulting  from  the  initial  angle 
of  attack,  and  from  the  first  and  second  twist  distributions  can  be  calcu¬ 
lated  from  the  following  equations 


R 


a 

■l 


v=i 


v- 


(B7) 
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3 

°!A  *  X  (ClXav 

v=i 


%  *  I 

v=l  v 

The  rolling-moment  coefficient  for  the  flexible  wing  for  any  desired 
angle-of -attack  distribution,  av,  is  then  obtained  by  substituting  values 
from  the  above  equations  into  equation  (B2).  The  reference  rolling -moment 
coefficients  with  the  aerodynamic  loads  at  the  quarter  chord  were  calcu¬ 
lated  through  use  of  the  previous  equations  (with  q  =  1  psi)  and  are 
tabulated  in  the  following  table: 


(B8) 


(B9) 


V 

1 

2 

3 

C 7 

LV 

0.120 

0.298 

0.326 

c, 

?AV 

-.0376 

-.0698 

-.0376 

Cl 

\ 

.OO856 

.01553 

.00796 

In  order  to  determine  the  rolling-moment  coefficient  for  a  specific 
derivative,  the  wing  angle-of -attack  distribution  must  be  known.  For  a 
unit  sideslip  angle,  the  wing  contribution  to  may  be  separated  into 

loadings  due  to  two  angle-of-attack  distributions  (ref.  13):  (l)  a  con¬ 

stant  spanwise  angle-of-attack  distribution  due  to  the  difference  in 
effective  velocity  acting  on  each  wing  panel  for  the  yawed  attitude  and 
(2)  an  angle-of-attack  distribution  proportional  to  the  wing  dihedral 
angle,  which  for  this  case  may  be  considered  to  be  due  entirely  to  the 
upward  deflection  of  the  wings  resulting  from  the  level  flight  symmetrical 
air  loads. 

For  the  rigid  airplane,  the  contribution  to  C,  due  to  the  differ- 

p 

ential  velocities  acting  on  each  wing  panel  was  estimated  from  reference  l4 

as  Cln  /CL  =  -°-129.  Substituting  values  for  a  unit  angle  of  attack  into 
Pi 

equations  (B7),  (B8) ,  (B9) ,  and  (B2)  gives 


30 


NACA  TN  3874 


-0.129  1  " 


0.14-51 

0.744 


0.0321  n 
0.1451  q‘ 


=  -0.129 


0.195Q  ^ 

1  +  0.221 qJ 


The  dihedral  angle,  measured  in  a  plane  parallel  with  the  YZ  plane, 
was  calculated  through  use  of  the  influence  coefficients  for  symmetrical 
loads  described  in  reference  1  and,  for  the  flight  range  of  interest,  could 
be  expressed  as 


TJ=1 


c 


L 


0.10351 


The  calculated  distribution  of  dihedral  angle,  normalized  with  respect  to 
the  value  at  the  tip,  is  given  in  the  following  table: 


4 

0.383 

0.707 

0.924 

r 

.520 

.900 

.992 

For  the  wing  at  a  unit  angle  of  sideslip,  these  values  of  dihedral  angle 
represent  changes  in  the  wing  angle  of  attack,  and  the  rolling-moment 
coefficient  due  to  dihedral  angle  may  be  determined  from  equations  given 
previously. 


-(0.10359) (0.557) 


0.1197 

0557 


i 


0.02664 

1  + - 

0.1197 


=  -0.05771 


0.2l4q  \ 
1  +  0 . 222q/ 


Since  only  first-order  effects  have  been  considered,  the  changes  in 
dihedral  angle  due  to  antisymmetric  loads  have  been  neglected. 
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The  total  value  for  the  wing  contribution  to  C,  becomes 


=  -0.129  1  - 


0.1951 
1  +  0.221q. 


-  0.05771  1  - 


0.2l4q 
1  +  0.222q.. 


(BIO) 


While  the  two  load  distributions  yield  somewhat  different  values  of  the 
ratio  C^/C^-r*  ^h e  values  of  k  are  essentially  the  same.  This  was 

also  found  to  he  true  for  other  types  of  load  distributions.  Note  that 
the  values  of  obtained  from  equation  (BIO)  are  not  the  same  as 

those  given  in  table  III  since  only  the  contribution  of  the  wing  has 
been  considered  here. 
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TABLE  I.-  PHYSICAL  CHARACTERISTICS  OF  THE  TEST  AIRPLANE 


Wing 

Area,  sq  ft  . 

Span,  ft . • . 

Aspect  ratio  .  .  . 

Taper  ratio  . 

Mean  aerodynamic  chord,  ft  . 

Sweepback  of  quarter- chord  line,  deg  . 

Airfoil  thickness  ratio  (parallel  to  body  center  line), 

percent  . 

Dihedral ,  deg . . 

Ailerons 

Area  aft  of  hinge  center  line  (each),  sq  ft  . 

Aileron  span  to  wing  span  ratio  . 

Average  aileron  chord  to  wing  chord  ratio  . 

Vertical  tail 

Area  (including  dorsal),  sq  ft  . 

Span,  ft  . 

Aspect  ratio  . 

Taper  ratio  ..  . 

Mean  aerodynamic  chord,  ft  . 

Sweepback  of  quarter- chord  line ,  deg  . 

Airfoil  thickness  ratio  (parallel  to  body  center  line), 

percent . . . 

Distance  l/4-chord  wing  M.A.C.  to  l/4-chord  vertical  tail 

M.A.C.,  ft  . 

Rudder 

Area  aft  of  hinge  center  line,  sq  ft  . 

Average  rudder  chord  to  vertical  tail  chord  ratio  .... 

Average  gross  weight,  lb  . 

Moment  of  inertia  about  the  X  axis ,  slug-ft2  . 

Moment  of  inertia  about  the  Z  axis,  slug-ft2  . 

Inclination  of  principal  longitudinal  axis  with  respect 
to  fuselage  reference  line,  deg  . 


1428 

116 

9.^3 

0.42 

13.0 

35 


53.8 

0.405 

0.26 

227 

18.7 

1.5^ 

0.3^ 

13.02 

3^-9 


51.2 

0.30 

115,000 

074,000 

306,000 


TABLE  II.-  STRUCTURAL  I  INFLUENCE  COEFFICIENTS 


front  spar  location;  R  denotes  rear  spar  location. 
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TABLE  III.-  LATERAL- DIRECTIONAL  STABILITY  DERIVATIVES  PREDICTED  FROM  WIND- 
TUNNEL  DATA  AND  AVAILABLE  THEORY;  h  =  35,000  FEET,  W  =  115,000  POUNDS 


Quantity 

Flexible 

Rigid 

M 

0.5 

0.6 

0.7 

0.8 

0.5 

0.6 

0.7 

0.8 

a 

87.2 

125.5 

171.0 

223.1 

87.2 

125.5 

171.0 

223.1 

CL 

.922 

.642 

.471 

.361 

.922 

.642 

.471 

.361 

% 

-.430 

-.429 

-.427 

-.426 

-.490 

-.507 

-.540 

-.570 

Cn 

nP 

-.0960 

-.0660 

-.0470 

-.0350 

-.0994 

-.0692 

-.0507 

-.0389 

c7 

ir 

.215 

.136 

.090 

.060 

.264 

.184 

.135 

.104 

cn 

ur 

-.14-2 

-.137 

-.134 

-.130 

-.150 

-.150 

-.150 

-.150 

% 

-.170 

-.141 

-.122 

-.112 

-.157 

-.129 

-.113 

-.101 

% 

.110 

.106 

.102 

.098 

.120 

.120 

.120 

.120 

% 

-.575 

-.566 

-.555 

-.544 

-.600 

-.600 

mm 

-.600 

C7 

Z*r 

.0052 

.0090 

.0109 

.0122 

.0057 

.0095 

.0119 

.0135 

Cnfi 

6r 

-.0724 

-.0700 

-.0672 

-.0645 

-.0778 

-.0778 

-.0778 

-.0778 

Cy 

Y5r 

.183 

.180 

.175 

.171 

.194 

.194 

.194 

.194 

C7 

*6a 

.0594 

.0528 

.0462 

.0392 

.0770 

.0770 

.0770 

.0770 

Cn„ 

5a 

-.0103 

-.0064 

-.0041 

-.0027 

-.0133 

-.0093 

-.0068 

-.0052 

rx 

-.116 

-.059 

-.005 

.031 

-.103 

-.047 

-.010 

.014 

rz 

-.048 

-.021 

-.005 

.007 

-.047 

-.020 

-.005 

.007 

38 


NACA  TN  387^ 


Figure  2.-  Two- view  drawing  of  test  airplane. 


ko 


NACA  TN  387^ 


t.  sec 


(t>)  Time  histories  of  (3,  rcg  ,  rt  ,  n^,  and 

B  h 

Figure  3.-  Concluded. 
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- - Experimental  data 


- Results  from  analog  computer 


t,  sec 


(a)  Time  histories  of  5r,  p  ,  p.  ,  and  |3. 

CSB  XB 

Figure  A.-  Typical  time  histories  for  rudder  pulse;  M 
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Figure  5.-  Wing  stations  and  weight  distribution  used  in  analysis 
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Figure  6.-  Influence  of  flexibility  on  effective  inertia  characteristics. 


Phase  angle 


Phase  angle,  <Pr  deg  Amplitude  ratio,  — per  sec 

cq  b  7 
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Frequency,  aj,  radians/sec 


(h)  Yawing  velocity  at  the  center  of  gravity. 
Figure  8.-  Continued. 


Phase  angle,  f  deg  Amplitude  ratio, 


MCA  TN  3874 


Phase  angle,  <£>  ,  deg  Amplitude  ratio,  — r-  ,  g/deg 
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Frequency,  oj,  radians/sec 
(d)  Wing- tip  acceleration 
Figure  8.-  Concluded, 


rGg  g 

Phase  angle,  <PD  ,  deg  Amplitude  ratio,  — ; —  ,  per  sec 
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Figure  9*“  Frequene 


Phase  angle,  <i>r  ,  deg  Amplitude  ratio,  9b  ,  per  sec 
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Figure  9.-  Continued. 


Phase  angle,  ,  deg  Amplitude  ratio,  per  sec 


MCA  TN  3874 
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Frequency,  w,  radians/sec 


(c)  Yawing  velocity  at  the  tail. 


t 


Figure  9.-  Continued. 


a 


Frequency,  cu,  radians/sec 
(b)  Yawing  velocity. 


Figure  11.-  Concluded 


Experiment 
Predicted  flexible 
Predicted  rigid 
Results  from  analog 
computer  matching 


Phase  angle,  <I>r  ,  deg  Amplitude  ratio,  — cq  B  per  sec 
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Frequency,  oj,  radians/sec 

(h)  Yawing  velocity  at  the  center  of  gravity. 


Figure  12.-  Continued 


Phase  angle,  <t>Tt  ,  deg  Amplitude  ratio,  |“H,  per  sec 
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O  Experiment 

-  Predicted  flexible 

- Predicted  rigid 
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Frequency,  cu,  radians/sec 


(c)  Yawing  velocity  at  the  tail. 


Figure  12.-  Continued 


Phase  angle,  ,  deg  Amplitude  ratio, 
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(d)  Sideslip  angle. 
Figure  12.-  Continued. 


Phase  angle,  <Pn  ,  deg  Amplitude  ratio,  — ,  g/deg 


MCA  TN  387^ 


Frequency,  cu,  radians/sec 
(e)  Wing-tip  acceleration. 


Figure  12.-  Concluded. 


Phase  angle,  D  ,  deg  Amplitude  ratio,  -  ,  per  sec 
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Frequency,  oj,  radians/sec 

(a)  Rolling  velocity  at  the  center  of  gravity. 

Figure  13.-  Comparison  of  experimental  and  predicted  frequency  responses 

to  rudder  input;  M  =  0.71. 


Phase  angle,  <Pr  ,  deg  Amplitude  ratio,  8  ,  per  sec 


MCA  TN  387^ 


Frequency,  a>,  radians/sec 

(c)  Yawing  velocity  at  the  center  of  gravity 
Figure  13.-  Continued. 
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(e)  Sideslip  angle. 
Figure  13.-  Continued. 
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amplitude  ratios 


Experiment 
Predicted  flexible 
Predicted  rigid 
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Figure  *15.-  Concluded. 
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NACA  -  Langley  Field,  Va. 


NACATN  3874  1-  Stability,  Lateral  -  NACA  TN  3874  1-  Stability, 

National  Advisory  Committee  for  Aeronautics.  Static  (1.8. 1.1. 2)  National  Advisory  Committee  for  Aeronautics.  Static 
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